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ABSTRACT 


Studies  have  continued  of  the  ability  of  a  single  Large  Aperture  Seismic  Array 
(LASA)  station  to  detect  and  roughly  locate  teleseismic  events  and  to  generate 
outputs  derived  from  various  long-  and  short-period  seismogram  parameters, 
particularly  those  related  to  exact  epicenter  location  and  source  discrimina¬ 
tion.  The  behavior  of  networks  that  exchange  such  data  and  that  consist  of 
both  small  and  large  array  stations  has  begun  to  receive  attention.  Computer- 
controlled  automatic  monitoring  and  maintenance  of  large  seismometer  arrays 
are  discussed. 
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INTRODUCTION 


This  is  the  sixth  Semiannual  Technical  Summary  Report  on  Lincoln  Laboratory's 
work  for  the  Advanced  Research  Projects  Agency  on  the  seismic  discrimination 
problem  (Vela  Uniform). 

Beginning  in  this  reporting  period,  we  have  attempted  to  shift  the  main  focus  of 
our  work  away  from  the  Large  Aperture  Seismic  Array  (LASA)  as  an  end  in  it¬ 
self  toward  consideration  of  the  seismic  surveillance  problem  as  a  whole  and  the 
role  that  networks  of  arrays  of  various  kinds  might  play  in  this  function. 

Studies  of  networks  (Sec.  I)  are  presently  limited  to  mathematical  simulations  and 
to  studies  of  LASA  data  as  sample  outputs  from  a  station  to  a  net  of  stations.  The 
particularly  interesting  outputs  are  various  waveform  discriminants,  which  are 
treated  in  Sec.  II. 

The  on-site  detection  and  location  that  take  place  autonomously  at  a  single  large 
array  station  are  described  in  Sec.  III.  These  detections  and  locations,  performed 
at  one  such  station,  can  be  used  to  initiate  off-line  array  processing  functions  at 
this  and  other  stations.  Section  IV  describes  recent  results  in  the  off-line  proc¬ 
essing  area,  both  with  short-  and  long-period  array  structures. 

The  current  program  on  techniques  of  automatic  maintenance  and  monitoring  of 
all  physical  elements  of  large  arrays,  along  with  recent  additions  and  improve¬ 
ments  made  to  the  Montana  LASA,  is  described  in  Sec.  V. 


P.  E.  Green 
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SEISMIC  DISCRIMINATION 


I.  SEISMIC  SURVEILLANCE  NETWORKS 


A.  LARGE  ARRAYS  IN  SEISMIC  SURVEILLANCE  NETWORKS 


During  this  reporting  period,  some  thought  was  given  to  the  role  that  one  or  more  large 

arrays  such  as  the  LASA  in  Montana  might  play  as  components  of  networks  to  detect,  locate  pmd 

identify  seismic  events.  Previous  experience  with  the  single  Montana  LASA  was  used  to  define 

the  outputs  that  large  and  small  array  stations  can  contribute  to  a  joint  network.  The  operation 

of  the  network  had  to  be  inferred,  partly  through  the  use  of  simulation  of  the  type  described  in 

Sec.  B  below.  u  ur  „  . 

H.  W.  Briscoe 

P.  E.  Green 

E.  J.  Kelly 


B.  SYSTEM  SIMULATION  OF  EPICENTER  LOCATION 

A  system  simulation  program  has  been  written  to  permit  comparison  of  various  world -wide 

networks  from  the  point  of  view  of  their  ability  to  provide  accurate  epicenters.  A  set  of  station 

coordinates  and  noise  levels  is  assigned,  and  it  is  assumed  that  epicenters  are  determined  from 

reported  arrival  times  by  a  least-squares  procedure.  This  procedure  weights  each  station  re- 

2 

port  inversely  with  the  variance,  c 7  ,  of  the  error  in  its  time-picking  (i.e.,  a  maximum-likelihood 
estimate  for  Gaussian  errors).  This  variance,  in  turn,  is  determined  from  the  formula 


or  =  a(A/  n)”1 


+  b 


where  (A/n)  is  the  voltage  signal -to-noise  ratio  at  the  station,  and  a  and  b  are  constants.  Jf 
A/ n  <  y,  or  if  the  station  is  beyond  100°  from  the  epicenter,  no  detection  is  assumed.  We  have 
used  the  values  a  =  1.0  and  b  =  0.1  (in  seconds)  so  that  time -accuracy  standard  deviations  range 
from  0.1  to  just  over  2.0  sec.  So  far,  we  have  treated  the  source  depth  as  known.  All  travel 
time  irregularities,  including  regional  source  biases  and  observing  station  corrections,  are  as¬ 
sumed  known  because  they  are  fixed  nonfluctuating  and  presumably  measurable  properties  of 
the  earth. 

A  linearized  error  analysis  of  this  estimation  procedure,  due  originally  to  M.  J.  Levin,  is 
used  to  obtain  the  parameters  of  the  error  ellipses  for  the  resulting  epicenters.  Because  station 
corrections  are  assumed  known,  the  estimates  are  unbiased.  Although  this  feature  is  unrealistic, 
the  intent  is  to  compare  various  networks  on  the  basis  of  the  random  component  of  their  errors 
of  epicenter  determination.  For  each  epicenter  on  the  preassigned  grid,  the  program  determines 
the  event  magnitude  required  so  that  the  order  of  the  ellipse,  of  95-percent  confidence,  is  an 
assigned  number  (we  have  used  5  00  km\  The  required  magnitude*  is  automatically  infinite  if 
fewer  than  three  stations  detect  the  event. 

Large  arrays  can  be  represented  in  such  a  study  in  various  ways.  For  one,  the  array  ban 
be  credited  with  a  certain  signal-to-noise  enhancement  factor,  y,  which  is  modeled  by  assigning 
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Section  I 


the  array  a  noise  level  smaller  by  the  factor  y  than  that  of  a  conventional  station.  Another  way 

is  to  enter  the  array  in  the  network  as  N  separate  stations,  where  N  is  the  number  of  subarrays 

2 

in  the  large  array.  This  automatically  reduces  a  for  the  array  by  the  factor  N,  which  is  a 

1/2 

reasonable  way  of  modeling  performance.  In  the  first  method,  y  is  taken  to  be  N  '  .  Results 

will  be  published  in  a  separate  report.  ^  j 

R.  J.  Kolker 
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II.  IDENTIFICATION 


A.  VISUAL  IDENTIFICATION  OF  DEPTH  PHASES 


In  an  experiment  reported  earlier,  Briscoe  and  Sheppard  found  that  the  "depth  phases"  pP 
and  sP  can  be  identified  with  high  reliability  at  an  individual  LASA  station.  Two  factors  other 
than  SNR  gain  appeared  to  contribute  (1)  space  diversity  whereby  the  P- waveforms  in  various 
subarrays  show  slight  differences,  particularly  in  the  coda,  but  have  the  depth  phase  in  common, 
and  (2)  the  array's  directivity,  the  array  aperture  being  large  enough  so  that  the  depth  phase 
shows  the  same  time  lag  after  P  from  subarray  to  subarray,  unlike  other  phases,  interfering 
teleseisms  or  reverberation. 

As  a  check  on  the  results  of  this  experiment,  it  was  decided  to  compare  the  record  in  cor¬ 
rectly  picking  depth  phases  achieved  by  the  Montana  array,  the  five  Vela  observatory  station$ 
individually  and  collectively,  and,  if  possible,  a  hypothetical  network  that  included  overseas 
single-sensor  stations.  (Each  Vela  observatory  is  roughly  equivalent  to  a  LASA  subarray.)  For 
this  purpose  we  picked  an  interval  of  time  (1  August  through  15  September  1966),  and  a  geograph¬ 
ical  region  (the  area  within  600 km  of  the  boundaries  of  the  Sino-Soviet  Bloc  area  and  within  90° 
of  LASA),  and  analyzed  the  64  events  so  reported  by  the  USCGS  ("Preliminary  Determination  of 
Epicenters")  from  world-wide  observations,  and  the  81  so  reported  in  the  LASA  station  bulletins. 
During  this  period  the  LASA  station  bulletins  were  prepared  from  film  records  of  nine  subarray 
outputs  and  no  beams  so  that  the  station's  performance  in  picking  depth  phases  was  well  below 
that  simulated  in  the  experiment  reported  earlier.  Figure  1  shows  the  results  of  the  analysis. 

Of  the  64  events  reported  by  CGS,  the  number  detected  by  each  observatory  is  given  as  the  over¬ 
all  height  of  a  corresponding  bar,  while  the  number  detected  by  at  least  one  observatory  in  a 
hypothetical  network  of  five  small  arrays  was  59  as  shown  by  the  next  bar.  LASA  detected  and 
located  43  of  the  64  events,  but  also  detected  and  located  (in  the  prescribed  region)  38  events 
not  reported  by  CGS;  this  total  of  81  is  indicated  by  the  bar  at  the  right. 

The  CGS-tabulated  station  reports  on  the  64  events  were  examined  to  see  how  often  each 
observatory  station  and  the  Montana  LASA  reported  a  phase  within  60  sec  after  P.  These  pP-P 
times  reported  were  then  checked  against  the  CGS-reported  depths.  If  they  lay  more  than  5  sec 
outside  the  CGS  one-sigma  tolerance  range  on  depth  computed  from  P-times  or  deduced  from 
depth  phases,  the  reports  were  called  false;  otherwise,  they  were  considered  to  be  valid  depth 
phase  identifications  (or  PcP  if  the  time  was  within  5  sec  of  theoretical  PcP  time).  The  number 
of  valid  identifications  and  false  reports  for  the  various  stations  are  shown  in  the  figure  as  shad¬ 
ing  within  the  vertical  bars.^ 


*H.  W.  Briscoe  and  R.  M.  Sheppard,  "A  Study  of  the  Capability  of  a  LASA  to  Aid  the  Identifica¬ 
tion  of  a  Seismic  Source,"  Technical  Note  1966-38,  Lincoln  Laboratory,  M.I.T.  (11  July  1966), 
DDC  437440,  H-731. 

tin  the  case  of  the  five-array  network,  a  false  report  was  counted  only  if  all  five  stations  re¬ 
ported  falsely.  In  the  case  of  LASA  it  was  assumed  that  the  38  events  on  which  no  CGS  report 
was  obtained  subdivided  into  false  ar*d  valid  reports  in  the  same  ratio  as  with  the  49  that  CGS 
did  report.  This  ratio  was  verified  in  an  independent  study  of  249  CGS-reported  (world-wide) 
events  on  which  LASA  or  one  or  more  observatories  reported  a  phase  within  60  sec  after  P. 
LASA  correctly  reported  157  depth  phases  or  PcP's  with  18  false  alarms;  the  five-array  net 
correctly  reported  125  depth  phases  or  PcP's  with  28  false  alarms. 


3 


Section  II 


Fig.  1.  Depth  phase  reporting  of  LASA  vs  a  5-station  network  in  the  United  States 
(1  August  -  15  September  1966;  distance  <  90°,  Bloc  countries  +6°). 
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Fig.  2.  Spectral  ratio  of  seismic  events. 
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Examination  of  CGS  reports  of  observations  of  these  same  events  from  single  sensor  seis 
mological  stations  around  the  world  revealed  that  no  one  such  station  was  nearly  so  effective 
either  detection  or  identification  of  depth  phases  as  just  one  of  the  Vela  observatories.  The 
single-sensor  stations  had  relatively  high  false  alarm  rates  on  depth  phase  reporting,  so  the 
difference  cannot  be  ascribed  to  insensitivity  of  record  analysis  at  the  station.  It  seems  cle&r 
that  the  greater  effectiveness  of  the  Vela  arrays  is  due  to  the  SNR  gain  inherent  in  the  multiple 
traces  plus  the  small  amount  of  space  diversity  available  within  their  3-  to  10-km  apertures 
Figure  1  shows  that,  even  using  only  nine  subarray  outputs  displayed  side  by  side  and  no 
directive  beams,  the  single  LASA  made  over  twice  as  many  valid  depth  phase  or  PcP  identifica¬ 
tions  as  the  five-station  network.  This  has  implications  in  organizing  systems  of  stations  and 

J 

deciding  whether  they  should  interchange  merely  station  bulletins  (lists  of  event  times)  or  wave¬ 
form  outputs.  Presumably,  the  large  disparity  in  depth  phase  performance  shown  in  Fig.  1 
would  not  exist  if  output  traces  from  the  five  observatories  had  been  transmitted  to  a  common 
point  and  displayed  side  by  side  for  analysis.  This  would  involve  problems  in  compensating  for 
propagation  delays  which  are  not  severe  for  real-time  work  with  the  200-km  LASA  aperture, 
but  which  would  become  formidable  for  a  five-array  net  if  many  epicenter  regions  were  to  be 
examined.  The  data  of  Fig.  1  corroborate  in  a  qualitative  way  the  encouraging  picture  of  the 
depth  phase  identifying  ability  of  a  single  LASA  given  by  the  earlier  results. 

P.  E.  Green 


B.  RATIOS  OF  SPECTRAL  DENSITIES 

The  possible  use  of  spectral  differences  between  explosive  and  natural  events  as  a  means 
of  identifying  the  nature  of  a  seismic  source  has  been  suggested  and  investigated  by  several 
people.  Near-zone  measurements  have  shown  that  explosions  generate  seismic  signals  rich  in 
high-frequency  energy  (about  1.5  Hz),  but  most  of  this  energy  is  filtered  out  in  propagation  to 
teleseismic  distances.  Most  investigations  of  the  short-period  P-energy  have  been  limited  by 
low  SNR  except  for  a  very  narrow  band  around  1  Hz. 

Since  one  of  the  potential  advantages  of  a  large  array  is  to  provide  a  good  SNR  on  signalj 
energy  over  a  wider  bandwidth  than  that  available  from  a  smaller  station,  a  study  is  in  progress 
to  investigate  whether  the  wider  bandwidth  will  help  detect  spectral  differences  between  explosive 
and  natural  sources.  The  initial  procedure  has  been  to  use  beamforming  (delay-and-sum)  proc¬ 
essing  to  improve  the  SNR  on  a  signal  and  then  to  compute  the  ratio  of  the  energy  in  a  band  from 
0.3  to  0.7  Hz  to  the  energy  in  the  band  from  1.5  to  1.9  Hz.  This  ratio  should  be  higher  for  earth¬ 
quakes  than  for  explosions.  A  set  of  27  events  including  explosions  from  four  underground  test 
sites  was  processed  using  this  procedure.  Results  of  this  preliminary  experiment,  which  are 
shown  in  Fig.  2,  tend  to  show  some  separation  of  the  two  source  types.  One  exception  is  apparent 
in  this  figure.  A  presumably  explosive  event  of  LASA  magnitude  4.9  has  a  ratio  of  0.44.  This 
was  found  to  be  the  result  of  an  unusually  high  level  of  low-frequency  noise  after  processing  so 
that  the  low-frequency  window  was  actually  looking  at  noise,  not  signal.  The  delay-and-sum 
processing  was  not  effective  in  this  case,  partly  because  data  from  several  subarrays  contained 
transients  due  to  phone-line  trouble  and  calibrations  being  performed  during  the  time  the  event 
arrived,  so  they  could  not  be  used  in  the  beamforming  process.  More  elaborate  processing 
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c 

Fig. 3.  Percent  of  earthquakes  with  complexity  ^ C . 
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techniques  are  being  investigated  to  provide  better  signal-to-noise  level,  particularly  when  onty 
part  of  the  array  is  operating  normally. 

The  small  sample  used  in  this  preliminary  experiment  was  chosen  as  a  pilot  population,  ahd 
a  much  larger  sample  including  all  events  from  the  experiment  described  in  Sec.  III-D  will  be 
processed  when  the  procedure  has  been  evaluated  and  debugged. 

H.  W.  Briscoe 

C.  COMPLEXITY  MEASUREMENTS 

The  automatic  station  bulletin  program  (discussed  in  Sec.  Ill -A)  computes,  among  other 
things,  the  complexity  of  each  detected  P-arrival  at  each  subarray,  and  averages  those  com¬ 
plexities  which  correspond  to  subarray  output  channels  used  in  the  determination  of  epicenter. 

As  used  here,  the  complexity  of  a  waveform  f(t),  is  defined  by 

/53°  I  f(t)  |  dt 

f5  )f(t)  |  dt 

''o 

where  t  =  0  represents  the  time  of  automatic  detection  (t  in  seconds).  For  a  single  event,  the 
measured  complexities  are  observed  to  vary  widely  across  the  array,  and  it  was  thought  that  the 
average  value  might  be  a  more  reliable  discriminant  than  single  trace  values.  Another  approach 
to  averaging  over  the  array  was  described  in  Sec.  V-B  of  the  last  Semiannual  Technical  Summary 
Report.* 

The  automatic  bulletin  program,  using  a  notch  prefilter  to  minimize  microseismic  noise, 
was  run  on  38  earthquakes  and  12  explosions  in  our  data  library.  The  earthquakes  were  all 
P-arrivals  (i.e.,  distance  not  greater  than  100°)  and  included  none  whose  reported  depth  exceeded 
200  km,  on  the  grounds  that  deeper  events,  while  often  simple  in  waveform,  would  be  recognized 
as  deep  from  P-arrival  times  alone.  The  explosions  were  from  three  source  points.  The  results 
are  shown  in  Fig.  3,  a  plot  of  percentages  of  the  two  populations  whose  complexities  exceeded  the 
value  of  the  abscissa.  None  of  the  explosions  in  this  group  had  a  complexity  greater  than  2.5, 
while  87  percent  of  the  earthquakes  were  more  complex  than  this.  Three  of  the  explosions  had 
complexities  in  excess  of  1.6,  and  these  were  the  three  weakest  explosion  records  in  the  group 
(the  largest  amplitude  was  limp).  Noise  alone  has  a  complexity  of  about  6,  and  weak  signals 
nearly  always  have  higher  complexities  than  they  would  have  had  at  high  signal-to-noise  ratios. 

All  of  the  earthquakes  in  our  population  had  complexities  in  excess  of  1.  6,  hence  the  power  of 
complexity  as  a  discriminant  could  perhaps  be  significantly  improved  by  computing  complexities 
only  on  array-processed  traces  of  high  SNR.  A  test  of  this  notion  is  in  progress. 

The  explosion  complexities  also  showed  a  relatively  small  variation  across  the  array.  Fop 
each  explosion  we  computed  the  standard  deviation  of  the  21  measured  complexities,  and  found 
that  the  average  value  (for  all  12  explosions)  was  0.4,  just  27  percent  of  the  average  complexity 
of  1.5.  The  earthquakes  showed  a  much  higher  variability. 


*  Semiannual  Technical  Summary  Report  to  the  Advanced  Research  Projects  Agency  on  Seismic 
Discrimination,  Lincoln  Laboratory,  M.I.T.  (30  June  1966),  DDC  637308. 
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The  recent  (27  October  1966)  event  in  Novaya  Zemlya  showed  strongly  anomalous  behavior 
for  a  presumed  explosion  (similar  to  former  events  from  that  region).  Its  complexity,  averaged 
over  the  array,  was  6.8,  with  a  standard  deviation  of  41  percent  of  that  value.  No  explanation  for 
this  is  available,  and  complexity  thresholds  may  ultimately  have  to  be  dependent  upon  epicenter. 

E.  J.  Kelly 


D.  RECTILINEAR  MOTION  STUDIES 

The  three-component  processor  described  in  Sec.  VII  of  the  last  Semiannual  Technical 
Summary  Report  has  been  operated  with  on-line  signals  telemetered  from  a  set  of  short-period 
seismometers  at  TFO  in  Arizona.  The  processor  is  basically  an  off-line  device  in  that  its 
parameters  (gains,  filter  center  frequencies  and  bandwidths)  ought  to  be  adjusted  and  optimized 
for  each  event  under  study.  However,  due  to  personnel  limitations  we  have  used  only  unattended 
real-time  film  monitoring  of  the  Arizona  data  with  parameters  set  to  form  the  best  compromise 
over  the  range  of  events  which  might  occur. 

One  expects  the  processor  to  enhance  an  analyst’s  ability  to  discover  and  identify  later 
phases  of  an  event,  including  pP,  out  of  the  large  amount  of  delayed  return  energy  arriving  after 
the  first  P.  The  processor  offers  an  advantage  over  operations  performed  on  exclusively  verti¬ 
cal  waveforms  in  that  it  is  sensitive  both  to  the  degree  of  rectilinearity  of  the  earth  particle  mo¬ 
tion  and  to  the  azimuthal  angle  of  arrival.  (Experience  has  shown  that  during  an  extended  coda 
of,  for  instance,  20-sec  duration,  75  percent  of  the  time  can  be  eliminated  as  possible  pP.  The 
remaining  5  sec  are  typically  split  up  among  several  discrete  peaks  which  the  processor  output 
suggests  are  possible  phases.)  It  is  a  difficult  problem  to  measure  quantitatively  the  amount  of 
improvement  in  phase-picking  the  processor  really  affords.  This  is  because  the  processor  is 
not  a  threshold  device  (such  as  an  event  detector)  and  because,  as  an  analyst's  aid,  its  output 
should  be  used  in  conjunction  with  other  analytical  techniques  (such  as  beam  outputs  steered  at 
different  velocities)  to  narrow  down  the  range  of  possibilities  in  making  actual  phase  picks. 

Nevertheless,  an  attempt  was  made  to  see  what  kind  of  phase-picking  was  possible  using  the 
output  of  the  present  processor  configuration  as  the  primary  indicator.  We  studied  all  available 
events  reported  by  the  CGS  during  the  period  1  August  through  15  October  1966  in  the  region 
within  6°  of  the  Sino-Soviet  Bloc  boundaries.  A  number  of  events  were  not  available  on  film; 
others  were  available,  but  of  too  large  or  small  an  amplitude  to  be  of  use  to  the  processor  with 
the  gain  settings  then  in  effect;  and  a  few  were  obscured  by  severe  man-made  noise  at  the  site. 
No  events  were  omitted  from  consideration  because  of  poor  waveform  or  difficulty  in  reading 
later  phases  where  the  initial  P-phase  was  visible. 

From  the  24  usable  events  so  obtained  (including  18  in  the  Kurile-Kamchatka  region),  an 
attempt  was  made  on  each  event  to  pick  the  single  cleanest  depth  phase  candidate  (except  once 
when  none  seemed  likely)  using  the  filmed  processor  output.  Two  general  principles  were 
followed:  (l)  a  depth  phase  should  present  a  processor  output  very  similar  in  appearance  to 
that  obtained  from  P  itself;  (2)  a  depth  phase  should  have  one  of  the  largest  subsequent  ampli¬ 
tudes  at  the  proper  azimuth.  These  two  criteria  for  phase-picking  occasionally  conflicted.  In 
such  cases,  it  was  usually  arbitrary  which  criterion  was  actually  applied.  The  time  of  the  phase 
candidate  so  selected  was  then  compared  with  possible  times  predicted  for  pP,  sP,  or  PcP  from 
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TABLE  1 

BLOC  AND  KURILE-KAMCHATKA  DEPTH  PHASE  IDENTIFICATIONS 

Bloc 

Kurile- 

Kamchatka 

No.  of  three-component  film  records  available 

37 

29 

No.  of  usable  three-component  records 

24 

18 

No.  of  which  apparent  depth  phase  selected 

23 

17 

No.  of  PcP's 

1 

1 

No.  of  false  alarms 

5 

3 

No.  of  pP  or  sP 

17 

13 

No.  of  usable  three-component  events  on  which 

TFO  reported 

24 

18 

No.  of  TFO-reported  depth  phases 

5 

5 

No.  of  PcP's 

0 

0 

No.  of  false  alarms 

0 

0 

No.  of  pP  or  sP 

5 

5 

the  CGS  location  and  depth  determinations  on  this  event.  If  the  selected  candidate  lay  within  the 
tolerance  range  for  pP  or  sP  (plus  or  minus  an  additional  5  sec  to  allow  for  timing  errors  in 
reading  the  records),  the  phase  pick  was  considered  successful.  Valid  PcP  picks  are  listed 
separately.  (These  are  the  same  criteria  used  in  the  experiment  described  in  Sec.  II-A.) 

Table  I  compares  our  TFO  three-component  readings  with  the  reports  to  CGS  from  TFO 
vertical  outputs  alone.  (Comparable  data  on  all  five  observatory  arrays  and  the  LASA  are  con^ 
tained  in  Fig.  1.)  It  is  seen  that  a  very  high  score  is  achieved  by  the  three-component  processor, 
accompanied,  however,  by  a  high  false-alarm  rate  of  about  one  to  every  four  valid  picks.  This 
is  to  be  contrasted  with  the  five-station  Vela  array  network  and  LASA,  which  achieved  a  com¬ 
parable  fraction  of  successful  picks  but  no  more  than  two  false  alarms  out  of  a  set  of  64  events* 
In  comparison  with  TFO  alone,  the  processor  gives  a  much  higher  rate  of  valid  picks  (three 
times  as  many)  but,  again,  a  much  higher  false  alarm  rate.  Although  it  cannot  be  shown  in 
Table  I,  which  is  based  on  a  single  pick  from  each  event,  the  processor  frequently  disclosed 
several  good  phase  candidates  and  it  frequently  occurred  that  the  second  and/or  third  candidates 
turned  out  to  be  valid  whether  or  not  the  first  candidate  itself  was  or  was  not  a  valid  phase. 

An  improved  three-component  processor  having  outputs  normalized  by  the  envelope  of  the 
vertical  seismogram  will  be  employed  in  further  studies. 

C.  A.  Wagner 
E.  Gehrels 
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m.  LASA  DETECTION  AND  LOCATION 


A.  LASA  DAILY  STATION  BULLETIN 

The  format  and  procedures  for  generating  the  LASA  station  bulletin  have  been  modified 
several  times  during  this  period.  In  August  the  format  was  modified  to  one  which  could  be  in¬ 
terpreted  automatically  and  was  easier  to  understand.  During  the  last  quarter,  an  automatic 
bulletin  generation  program  was  used  and  its  reliability  analyzed  to  suggest  areas  where  it  could 
be  improved.  In  November  additional  off-line  programs  to  aid  in  generation  of  the  bulletin  were 
put  into  routine  use  and  the  format  for  some  of  the  parameters  was  modified. 

Studies  of  the  amplitude  scatter  at  various  subarrays  in  the  array  are  summarized  in 
Table  II.  The  large  scatter  between  subarrays  compared  to  that  within  a  subarray  has  led  to 
the  conclusion  that  the  magnitude  estimate  for  an  event  obtained  by  averaging  over  all  of  LASA 
has  a  significantly  smaller  scatter  than  that  from  a  single  small  array  station.  As  a  result, 
one  of  the  changes  in  bulletin  format  in  November  was  to  report  an  average  magnitude  to  the 
nearest  0.1  magnitude  unit  instead  of  using  a  single  subarray  and  reporting  to  0.5  magnitude  uhit. 


TABLE  II 

COMPARISON  OF  P-AMPLITUDE  SCATTER 
OVER  7-  AND  200-KM  APERTURES 
(25  Events) 


Ratio  of  Standard  Deviation 
to  Mean 

Ratio  of  Strongest 
to  Weakest  Sensor 

Typical  Event 
(percent) 

Worst-Case 

Event 

(percent) 

Typical  Event 

Worst-Case 

Event 

Subarray  (16  sensors  in  F4) 

15 

30 

1.8:1 

2.6:1 

Array  (center  sensors 

from  each  of  16  subarrays) 

40 

80 

4:1 

9:1 

A  preliminary  evaluation  based  on  one  month  of  data  indicates  that  the  routine  use  of  the 
off-line  programs  designed  to  aid  the  preparation  of  the  bulletin  (including  the  initial  automati 
bulletin  program  described  in  Sec.B  below)  has  improved  the  threshold  for  event  reporting  in 
the  bulletin  by  0.2  to  0.3  magnitude  unit  compared  to  bulletin  preparation  without  machine  aid. 
During  normal  operating  conditions,  magnetic  tape  recordings  of  array  data  are  not  available 
for  24  hours  a  day,  7  days  a  week,  so  the  improved  threshold  will  only  be  realized  part  time. 

The  performance  of  the  bulletin  during  August  and  September  is  summarized  in  Table  III. 

The  total  improvement  over  the  last  6  months  in  the  overall  event  reporting  ability  of  the 
LASA  station  can  be  deduced  from  Fig.  4,  which  indicates  at  the  present  time  a  4.1  magnitude 
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TABLE  III 

LASA 

STATION  BULLETIN  ERRORS 

(Based  on  the  58  Events  Detected  and  Located  by  Both  LASA  and  CGS  as  Lying  ^  90° 
from  LASA  and  Within  6°  of  Sino-Soviet  Bloc;  1  August  —  15  September  1966) 

Magnitude  (LASA -CGS) 

Mean:  +0.09 

Standard  deviation:  0.31 

Epicenter  location 

Standard  deviation:  3.4°  (370  km) 

Origin  time 

Standard  deviation:  21.8  sec 

Fig.  4.  Number  of  events  detected  and  located  at  and  above 
a  given  magnitude  vs  that  magnitude. 
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threshold  above  which  75  percent  of  all  events  40°  to  90°  from  the  LASA  are  detected  and  located. 
The  dotted  curve  shows  results  obtained  in  the  spring  of  1966*  using  visual  operations  only, 
before  the  initiation  of  the  daily  bulletin. 

So  far,  the  bulletin  has  been  prepared  using  subarray  outputs  separately,  that  is,  without 

combining  them  into  beams  as  in  Secs.  Ill -C  and  -D.  Very  recently,  we  have  started  adding 

beam-detected  and  located  events.  TT  ,xr  _  . 

H.W.  Briscoe 


B.  AUTOMATIC  STATION  BULLETIN 

A  Fast  Automatic  STAtion  BULletin  (FASTABUL)  program  has  been  written  which  automat¬ 
ically  computes  and  tabulates: 

(1)  The  amplitude  of  the  first  arrival, 

(2)  The  period  of  the  first  arrival, 

(3)  The  complexity  of  the  seismic  waveform, 

(4)  The  direction  of  first  motion, 

(5)  The  horizontal  phase  velocity,  and 

(6)  The  azimuth  of  the  best -fitting  plane  wave  across  the  whole  array. 

From  these  last  two  quantities,  the  station  bulletin  also  computes: 

(7)  The  distance  of  the  epicenter  from  LASA, 

(8)  The  latitude  of  the  epicenter, 

(9)  The  longitude  of  the  epicenter, 

(10)  The  arrival  time  of  the  first  motion  at  LASA,  and 

(11)  The  origin  time . 

From  items  (1),  (2),  and  (7)  it  computes: 

(12)  The  magnitude  of  the  event. 

The  basic  problems  in  an  automatic  station  bulletin  are  to  have  the  computer  recognize  a 
P  (or  PKP)  first  arrival  from  the  background  noise  and  to  measure  accurately  these  arrival 
times  at  each  subarray.  This  program  goes  through  a  series  of  operations  on  the  seismic  wave¬ 
forms  and  picks  arrival  times,  culls  these  arrival  times  and  discards  the  ones  that  disagree 
with  a  plane  wave  propagating  across  the  array,  and  computes  the  station  bulletin  parameters. 

It  also  plots  the  seismograms  and  marks  on  them  the  times  it  has  picked.  Under  normal  con¬ 
ditions  it  does  rather  well,  but  sometimes,  as  when  the  signal  is  extremely  weak  or  when  two 
\  teleseisms  arrive  together,  it  does  poorly.  In  these  cases,  a  human  operator  can  look  over 
the  time  picks  apd  perhaps  make  a  wiser  choice  of  channels  for  the  computer  to  use  to  recom¬ 
pute  the  station  bulletin  parameters.  This  is  called  the  "semiautomatic  mode"  of  operation. 

The  input  waveforms  to  this  program  can  be  the  twenty-one  "10"  seismometer  outputs,  or 
21  maximum -likelihood  processed  traces  where  each  trace  is  a  combination  of  the  entire  jout- 
put  from  each  subarray,  or  waveforms  from  any  type  of  S/N  improvement  scheme.  Therefore, 
the  minimum  magnitude  for  which  this  automatic  station  bulletin  will  work  depends  on  the  type 
of  preprocessing  used  on  the  input  waveforms.  If  only  one  waveform  is  input  to  the  program 
(for  example,  a  beam  made  up  from  the  entire  array),  the  location  of  the  source  is  determined 


*  E.  J.  Kelly,  "LASA  On-Line  Detection,  Location  and  Signal-to-Noise  Enhancement,"  Technical 
Note  1966-36,  Lincoln  Laboratory,  M.I.T.  (1  July  1966),  DDC  636144,  H-727. 
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Fig.  5.  Error  in  location  (automatic 
station  bulletin). 


Fig.  6.  Error  in  magnitude. 


Fig.  7.  Error  in  origin  time. 
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by  where  the  beam  is  pointed,  so  the  program  need  only  calculate  complexity,  magnitude,  fii^st 
motion,  and  origin  time. 

The  exact  method  by  which  the  automatic  station  bulletin  parameters  are  computed  is  deh 
scribed  in  detail  in  a  recent  report.* 

This  program  was  tested  in  the  fully  automatic  mode  using  70  events,  which  are  all  the 
events  in  our  library  from  1  October  1965  to  13  January  1966  that  have  been  identified  with  CJGS 
events  and  for  which,  therefore,  we  know  the  location,  magnitude,  and  origin  time.  This 
should  represent  a  typical  population  of  events.  The  events  which  were  not  saved  in  our  library 
were  later  arriving  phases  of  strong  events. 

The  results  of  the  test  for  all  the  events  with  epicenters  less  than  110°  from  LASA  (P  phase 
only)  are  shown  in  Figs.  5,  6,  and  7.  The  error  is  plotted  as  a  function  of  LASA  amplitude  be¬ 
cause,  as  the  signal  becomes  stronger,  the  automatic  identification  becomes  better.  The  oyer- 
all  mean-square  error  using  fully  automatic  location,  averaged  over  events  of  all  sizes,  was 
6.2  great  circle  degrees.  The  mean-square  error  for  all  the  events  which  were  received  with 
an  amplitude  greater  than  2  mfi  was  4°.  These  can  be  reduced  somewhat  by  employing  the  semi¬ 
automatic  mode,  but  there  is  a  limit  of  about  2.5°  due  to  the  limited  aperture  of  the  LASA  airray 
itself  (200km)  and  inherent  timing  errors  which  include  our  present  imperfect  knowledge  of 
the  station  corrections.  These  timing  errors  become  much  more  restrictive  for  the  PKP  phases. 
The  errors  due  to  the  travel  time  curves  seem  small  compared  to  these  timing  errors,  evejn 
though  a  depth  of  33  km  is  assumed  throughout.  The  overall  mean-square  difference  between 
the  CGS  and  LASA  magnitudes  is  0.4,  but  this  is  not  too  surprising  since  the  CGS  reports  a[re 
averaging  stations  with  world -wide  coverage  and  which  have  large  deviations  themselves.  The 
overall  mean-square  error  in  determining  the  origin  time  is  33  sec. 

It  is  difficult  to  put  a  lower  bound  on  the  magnitude  for  which  this  automatic  station  bulletin 
will  work.  If  the  signal  is  received  with  an  amplitude  greater  than  5  mfji,  the  automatic  picjks 
are  made  about  as  well  as  a  human  can  make  them. 

Recently,  the  FASTABUL  program  has  been  modified  to  preprocess  the  individual  wave¬ 
forms  which  are  input  to  the  program.  The  preprocessing  used  here  means  passing  the  data 
\  through  a  four-pole  four -zero  filter  which  forms  two  notches  in  the  frequency  domain  at  0 jz 
'and  0.3  Hz.  The  filter  has  very  little  phase  shift  at  frequencies  higher  than  these  so  that  *  ng- 
Pjeriod  microseisms  are  attenuated  without  significantly  distorting  the  teleseism  waveforms. 

This  new  program  (NOTCHSTABUL)  has  been  tested  using  the  same  70  events  mentioned  above. 
/  Although  the  results  were  not  reduced  in  the  same  manner,  it  appears  that,  qualitatively,  the 
prefiltering  improves  the  performance  by  reducing  the  overall  mean-square  error. 

P.  L.  Fleck 


C.  OBSERVATIONS  IN  HONSHU  AREA 

The  difficulty  of  obtaining  a  direct  verification  of  the  detection  capability  of  LASA,  when 
some  form  of  predetection  processing  (such  as  beamforming)  is  employed,  owing  to  the  lack 


*  P.  L.  Fleck,  "FASTABUL  (A  Fast  Automatic  STAtion  BULletin  Program),"  Technical  Ndte 
1967-3,  Lincoln  Laboratory,  MILT.  (6  January  1967). 
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of  independent  reporting  of  the  requisite  weak  events,  has  frequently  been  pointed  out.  In  order 
to  obtain  reports  of  weak  events  in  time  to  permit  the  saving  of  tape  recordings  in  Montana,  we 
focused  our  on-line  beams  on  Japan  and  made  arrangements  with  the  Earthquake  Research 
Institute  of  the  University  of  Tokyo,  who  kindly  allowed  one  of  our  people  (J.  Fairborn,  a  grad¬ 
uate  student  in  the  M.I.T.  Geophysics  Department)  to  work  there  for  a  time  with  direct  access 
to  their  data.  The  data  itself  consisted  of  Helicorder-type  recordings  from  a  tripartite  of  sta¬ 
tions  surrounding  Tokyo.  These  stations  are  equipped  with  short-period  instruments  (peak 
response  at  20  Hz),  designed  for  the  detection  of  local  events.  The  times  and  amplitudes  of  the 
P  and  S  arrivals  for  a  5 -week  period  in  midsummer  were  recorded,  and  the  event  times  tele- 
typed  to  this  Laboratory.  As  a  result,  tape  recordings  of  about  150  events  were  obtained,  more 
than  half  of  which  are  from  the  Matsushiro  area. 

The  great  majority  of  these  events  were  within  10°  of  Tokyo,  and  we  developed  a  simple 
plane-earth  method  for  finding  the  epicenters  from  the  arrival  times  and  S-P  intervals.  The 
time-difference  between  the  arrivals  at  two  stations  implies  (using  an  empirical  P^-speed  of 
7.5km/sec)  that  the  epicenter  is  on  a  certain  hyperbola,  and  the  S-P  interval  for  one  of  these 
stations  fixes  a  distance,  yielding  two  possible  epicenters.  The  empirical  relation 


D  =  (7.86  +  0.04t )  r 


where  t  is  the  S-P  interval  in  seconds,  was  used  to  obtain  the  distance  D  in  kilometers.  If 
a  third  station  arrival  was  recorded,  the  ambiguity  between  the  two  epicenters  can  be  resolved. 

In  most  cases,  more  than  one  S-P  interval  was  recorded  and  often  all  three  arrival  times  were 
measured.  The  consistency  of  the  epicenters  obtained  from  such  redundant  data  was  fair,  and 
we  estimate  that  our  epicenters  are  accurate  to  ±15  km  for  most  of  the  events. 

The  next  step  is  the  off-line  formation  of  beams  steered  for  the  computed  epicenters  which, 
if  the  epicenters  themselves  are  to  be  believed,  will  establish  the  detection  threshold  for  beam¬ 
forming  in  the  optimum  case  where  the  events  occur  exactly  in  the  beam.  In  parallel  with  this, 
a  grid  of  14  beards  covering  Honshu  were  formed  off-line  on  all  recordings.  These  beams  are 
4°  apart  in  each  direction,  and  should  give  us  some  idea  of  how  future  beamsplitting  techniques 
will  work  out.  The  experiment  is  also  a  severe  test  of  our  station  corrections  for  this  area,  / 


and  should  provide  an  absolute  calibration  of  beam  pointing  of  considerable  precision. 


Although  signal  amplitudes  at  the  tripartite  stations  were  recorded,  we  do  not  have  suffi-/ 
cient  confidence  in  our  knowledge  of  the  effective  magnification  of  the  instruments  or  the  pre-  > 
cise  nature  of  the  first  arrivals  to  be  able  to  determine  (teleseismic)  body  wave  magnitudes. 
Hence,  amplitudes  at  LASA  are  being  measured  on  those  events  detected,  and  a  cumulative 
curve  of  number  of  detections  vs  apparent  magnitudes  is  being  constructed. 

The  on-line  beams  were  formed  using  rough  and  incomplete  station  corrections,  and  their 
performance,  as  read  from  film  recordings,  was  not  impressive.  However,  the  beams  formed 
off-line  have  used  station  corrections  for  the  entire  array  for  the  Honshu  area,  and  they  have 
shown  the  full  S/N  enhancement  factor  (*V"N)  that  has  been  found  with  beams  using  hand-picked 
times. 


*  "Travel-Time  Anomalies  at  LASA,"  Report  No.  LL-1,  Earth  Sciences,  A  Teledyne  Company, 
Alexandria,  Virginia  (2  November  1966). 
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TABLE  IV 

KURILES  UNDERWATER  TEST  EXPLOSIONS 

5  Ton 

1  Ton 

Total  number 

10 

7 

Data  gap  on  beams 

3 

3 

Visually  detected 

5 

1 

Possible  detection 

1 

1 

Not  detected 

1 

2 
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Section  III 


Reasonable  agreement  has  been  found  between  the  epicenters  determined  directly  from  the 
ERI  data  and  those  inferred  from  the  grid  of  beams.  Based  on  the  10  to  15  events  of  marginal 
visibility  we  have  so  far  found,  we  would  place  the  50-percent  probability  of  detection  threshold 
for  beams  at  about  magnitude  3.5  (fixed  level,  not  cumulative  probability). 

In  Fig.  8  we  show  an  example  of  the  Japan  data,  an  event  of  apparent  magnitude  3.7.  The 
upper  two  traces  are  straight  subarray  sums,  unfiltered,  and  the  next  four  are  beams,  also 
unfiltered,  steered  at  points  near  the  epicenter  (determined  from  ERI  data)  as  indicated  in  the 
figure.  These  beams  are  steered  sums  of  the  21  straight-sum  subarray  outputs.  We  are  also 
studying  the  performance  of  beams  made  up  of  seven  sensors  from  each  of  the  21  subarrays, 
each  sensor  appropriately  delayed,  and  using  various  sets  of  station  corrections. 

E.  J.  Kelly 
R. M.  Sheppard 


D.  KURILE  OBSERVATIONS 

During  November  and  December  1966,  as  part  of  another  Vela  Uniform  program,  a  network 
of  ocean-bottom  seismometers  was  placed  off  the  Kurile  Islands  and  a  series  of  17  high-explosive 
shots  were  fired  within  the  network.  In  order  to  take  advantage  of  the  local  coverage  of  small 
magnitude  events  provided  by  this  network,  all  five  beams  available  in  the  on-line  LASA  proc¬ 
essing  facility  were  steered  to  the  Kurile  region  and  an  attempt  was  made  to  operate  the  site 
24  hours  a  day.  Data  tapes  were  saved  for  all  events  detected  and  located  within  5°  of  the  Sino- 
Soviet  region  from  1  November  to  20  December,  and  for  all  events  detected  but  not  located  from 
15  November  to  20  December.  Events  using  the  beams  were  detected  by  visual  scanning  of  film 
records  of  the  beam  data.  Although  all  of  the  data  has  not  been  assembled  or  analyzed,  pre¬ 
liminary  tabulation  of  the  first  month  of  data  shows  44  located  Sino-Soviet  events  (within  5°  of 
the  Sino-Soviet  region)  between  1  November  and  6  December,  and  67  events  detected  but  not 
located  between  15  November  and  6  December.  Off-line  processing  of  several  days'  unlocated 
events  at  Lincoln  Laboratory  indicate  that  approximately  20  percent  of  the  currently  unlocated 
events  will  be  located  in  the  Sino-Soviet  region. 

The  results  of  visual  detection  of  the  17  high- explosive  shots  from  film  records  are  shown 
in  Table  IV.  The  film  records  were  unfiltered  or  filtered  with  a  wide-band  filter.  This  pres¬ 
entation  of  the  data  was  thus  not  designed  to  optimize  detection.  The  best  recorded  shot 
(3  December  at  22:49  GMT)  had  a  period  of  0.8  sec  and  a  LASA  magnitude  of  approximately  4.4. 

H.W.  Briscoe 
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IV.  ARRAY  PROCESSING 

A.  STUDIES  OF  OPTIMUM  SUBARRAY  GEOMETRIES  AND  PROCESSING 

1.  Implications  from  Beamforming 

In  an  earlier  report,*  some  results  were  described  which  showed  that  in  the  0.6-  to  %.0-Hz 
P-wave  spectral  range,  an  increase  of  minimum  inter-element  spacing  from  0.5  to  3  km  resulted 
in  a  significant  increase  in  SNR  gain  using  delay-and-sum  (DS)  processing  (beamforming).  The 
gains  achieved  were  roughly  the  same  as  those  available  from  maximum-likelihood  (FS)  process¬ 
ing.  These  results  argue  in  favor  of  redesigning  the  present  25-element  LASA  subarrayS  to  have 
an  aperture  of  around  20km  with  roughly  3-km  spacing  (see  Sec.V-C).  A  LASA  composed  of  such 
subarrays  would  develop  a  large  fraction  of  its  total  available  SNR  gain  by  on-line  beamforming. 

In  order  to  determine  an  optimum  configuration  for  the  design  of  future  LASA  subarrays,  a 
number  of  beamforming  experiments  were  carried  out,  using  LASA  data,  which  essentialiy  con¬ 
tinue  the  earlier  experiments.  The  LASA  data  were  prefiltered  with  the  special  0.6-  to  2.0-Hz 
Butterworth  filter  of  Fig.  9,  and  various  arrangements  of  seismometers  were  considered.  The 
results  are  presented  in  Table  V. 

It  is  seen  from  the  noise  reduction  figures  that  if  the  sensors  are  spaced  about  6  km  apart, 
the  noise  is  essentially  random  in  the  0.6-  to  2.0-Hz  frequency  range.  At  this  spacing,  the  noise 
reduction  is  13.2  db,  which  is  close  to  the  14  db  that  would  be  obtained  if  the  noise  were  inco¬ 
herent.  This  noise  reduction  is  consistent  with  the  noise  coherency  measurements  discussed 
in  the  next  section.  As  a  practical  matter,  it  is  inconvenient  to  have  6-km  spacing  between  sen¬ 
sors,  so  that  3-km  spacing  seems  more  reasonable.  If  25  sensors  are  spread  over  a  22-km  ap¬ 
erture  with  3-km  spacing,  an  average  SNR  gain  of  about  lldb  is  achieved.  If  this  minimum 
spacing  is  maintained  but  more  sensors  are  added  in  the  desire  to  increase  the  gain,  then  the 
aperture  of  the  subarray  must  be  increased.  However,  as  this  is  done,  the  signals  tend  to  be 
less  coherent  so  that  the  loss  in  gain  due  to  signal  reduction  increases. 

The  highest  average  SNR  gain  was  obtained  with  58  sensors  spread  over  40km.  However, 
the  loss  relative  to  n/"n  is  4.  7  db,  which  compared  unfavorably  with  the  3.0-db  loss  when  25  sen¬ 
sors  were  spread  over  22km.  Thus,  the  results  indicate  that  the  most  reasonable  choice  for 
the  design  of  a  future  LASA  subarray  is  to  have  25  sensors  with  a  minimum  3-km  spacing  be¬ 
tween  sensors  and  with  an  aperture  of  approximately  20  km. 

2.  Coherency  Measurements  on  Noise 

The  manner  in  which  the  coherency  of  seismic  noise  decreases  with  distance  provided  an 
indication  of  how  much  gain  can  be  obtained  by  delayed  sum  processing.  This  relationship  will 
be  developed  in  a  forthcoming  report. 

The  main  signal  energy  in  the  LASA  short-period  seismometer  is  at  frequencies  betWeen 
0.6  and  2.0  Hz.  Values  of  |C^(f)|  were  computed  at  five  frequencies  in  this  range  to  determine 


*  J.  Capon,  R.  J.  Greenfield  and  R.T.  Lacoss,  "Off-Line  Signal  Processing  Results  for  the  Large 
Aperture  Seismic  Array,"  Technical  Note  1966-37,  Lincoln  Laboratory,  M.I.T.  (11  July  1^966), 
PDC  637016,  H-730. 
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Fig.  9.  Gain  characteristics  of  various  filters. 
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TABLE  V 

RESULTS  OF  BEAMFORMING  EXPERIMENTS 

Aperture 

(km) 

Minimum 

Spacing 

(km) 

Number 

of 

Number 

of 

Noise 

Reduction* 

(db) 

Average 

Noise 

Reduction* 

Signal 

Reduction! 

(db) 

Average 

Signal 

Reduction! 

(db) 

Average 

SNR 

Ga  i  n 

/N 

Sensors  =  N 

Runs 

Min 

Max 

(db) 

Min 

Max 

(db) 

(db) 

7 

0.5 

25 

30 

4.0 

6.5 

6.1 

0 

0.2 

0.1 

6.0 

14 

22 

3 

25 

5 

10.0 

12.0 

11.5 

0.2 

0.7 

0.5 

11.0 

14 

40 

6 

25 

7 

10.9 

15.4 

13.2 

0.4 

2.5 

1.5 

11.7 

14 

40 

3 

58 

20 

11.6 

16.1 

14.5 

0.6 

4.2 

1.6 

12.9 

17.6 

65 

3 

91 

6 

14.3 

16.5 

15.3 

1.2 

6.1 

3.1 

12.2 

19.6 

*  Measured  with  respect  to  average  sensor  noise  power. 

t Measured  with  respect  to  average  sensor  signal  amplitude. 
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Fig.  10.  Noise  coherency  vs  seismometer  separation. 
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how  |C.j(f)J  varies  as  a  function  of  seismometer  separation.  The  value  of  |C.^(f)|  for  each  pair 
of  seismometers  in  a  LASA  subarray  was  computed  from  the  definition 

|C..(f)|  = 


l  fij(f) I 


U 


J f..(f)  f.(f) 
^  n  JJ 


where 


f  =  frequency 

f;,(f) 


C i j (f )  =  =  complex  coherency  for  i  and  seismometers 


f.j(f)  =  cross  power  spectra  for  i  and  j 


th 


seismometers. 


The  estimation  of  f_(f)  was  made  using  the  direct  segment  method.  The  data  block  length  used 
was  10  sec  which  gives  a  spectral  resolution  of  0.1  Hz.  The  length  of  data  used  for  estimation 


varied  from  410  to  700  sec  or  from  41  to  70  data  blocks. 


In  a  subarray,  the  distance  between  seismometers  varies  from  0.5  to  7  km.  This  range  of 
separations  was  divided  into  nine  subranges,  and  an  average  was  formed  for  each  subrange  using 
all  separation  pairs  in  the  subrange.  The  average  for  the  k^  subrange  is  denoted  by  C^(f).  The 
measurements  of  C^(f)  are  biased,  since  a  finite  length  of  record  was  used  for  the  spectral  es¬ 
timation.  The  bias  is  especially  large  for  values  of  C_(f)  near  zero.  A  correction  for  the  bias 
was  made  to  C^ff)  using  tables  given  by  Amos  and  Koopmans;^  the  number  of  degrees  of  freedom 
for  the  measurements,  called  n  by  Amos  and  Koopmans,  is  equal  to  the  number  of  data  blocks 
used  in  the  estimation.  The  corrected  values  of  C^(f)  are  denoted  by  C^ff).  In  Fig. 10,  CJ  (f)  is 
plotted  vs  the  separation  distance  covered  by  the  subrange.  The  results  show  that  the  noise  is 
more  coherent  at  the  lower  part  of  the  signal  frequency  band.  At  0.7  and  1.0  Hz,  C^(f)  is  gen¬ 
erally  below  0.2  for  separations  greater  than  3  km. 


3.  Convolutional  Filter  for  Short- Period  Data 

A  linear-phase  high-pass  convolutional  digital  filter  may  be  used  to  remove,  with  alrpost  no 
distortion,  microseismic  noise  while  passing  a  teleseismic  signal.  An  input  signal  X(t.)  ^s  fil¬ 
tered  by  the  following  formula  to  give  the  output  trace  Y (t . ) : 


2L 

Y(t.)  =  X(tj-  tQ)  -  Yj  X(t.  -  ktg)  •  A(ktg)  •  R(ktg) 
k=0 

where 


t  =  sampling  rate 
t 

L  =  =  half  the  length  of  the  filter 

s 


*  J.  Capon,  R.  J.  Greenfield  and  R.  J.  Kolker,  "A  Frequency-Domain  Synthesis  Procedure  for 
Multidimensional  Maximum-Likelihood  Processing  of  Seismic  Arrays,"  Technical  Note  1^66-29, 
Lincoln  Laboratory,  M.I.T.  (6  May  1966),  DDC  634233,  H-724. 

t  D.  E.  Amos  and  L.  H.  Koopmans,  "Tables  of  the  Distribution  of  the  Coefficient  of  Coherence  for 
Stationary  Bivariate  Gaussian  Processes,"  Monograph  SCR -483,  Sandia  Corporation  (March  1963). 
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RAW  OATA 


fc  =  0.5Hz,  121  POINT  FILTER 


fc=0.6Hzt121  POINT  FILTER 


fc*0.7Hz,121  POINT  FILTER 


fc  =  0.8  Hz,  121  POINT  FILTER 


fc  =  0.6  Hz  ,  201  POINT  FILTER 


Fig.  1 1 .  Example  of  convolutional  filter  of  Fig.  9  applied  to  record  of  20  November  1965  Marianas  event. 
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R(ktg)  =  sin 


»(kts  -  V 


A (kt  )  =  1  - 


kt  -  t 

1  o  o 1 


f  e  low-frequency  cutoff 

The  frequency-domain  response  of  this  filter  for  a  particular  choice  of  parameters,  as  well  as 
the  response  of  two  other  filters  used  in  LASA  data  processing,  is  shown  in  Fig.  9.  In  Fig.  ll, 
the  results  of  applying  the  convolutional  filter  to  seismic  data  show  how  little  this  filter  distorts 
a  teleseismic  signal. 

B.  LONG- PERIOD  RESULTS 

Preliminary  processing  of  data  from  the  21  three-component  long-period  seismometers  in 
LASA  has  been  completed.  (The  nominal  transfer  characteristic  of  the  seismometers  is  shown 
in  Fig.  18,  Sec.V.)  Delay-and-sum,  optimum  one-point  filters  (WDS)  and  maximum-likelihood 
filters  (FS)  have  been  applied  to  each  of  the  three  components  of  the  LP  data  from  LASA.  The 
processors  were  adjusted  to  pass  either  S  or  P  waves  from  teleseismic  events. 

We  have  reached  the  following  tentative  conclusions  based  upon  preliminary  results:  (1>  So¬ 
phisticated  linear  array  processors  can  be  applied  in  order  to  enhance  signals  on  long-period 
seismograms  with  no  measurable  signal  degradation  on  the  output  trace.  (2)  Voltage  SNR  gains 
in  the  signal  band  at  least  as  large  as  nTn  can  be  achieved  on  all  three  components  using  LP  (data, 
where  N  is  the  number  of  seismometers.  (3)  The  LP  noise  can  become  more  organized  causing 
SNR  gain  of  delayed  sum  processing  to  be  reduced  a  few  db  below  \T N,  but  not  adversely  affecting 
the  SNR  gain  of  optimal  processors.  (4)  Experiments  with  two  events  from  the  same  region  over 
a  week  apart  have  indicated  that  filters  which  are  designed  using  one  of  these  events  can  be  Used 
to  achieve  20-  to  25-db  signal  rejection  when  applied  to  the  other.  This  ability  to  see  through 
events  will  be  of  great  value  when  dealing  with  LP  data  because  of  the  long  duration  of  surface 
wave  trains. 

Table  VI  shows  typical  SNR  improvements  obtained  for  a  few  events.  The  three  forms  of 
processing  were  applied  to  three  strong  events  to  determine  SNR  gain  when  the  processor  i£  de¬ 
signed  to  pass  either  P  or  S  phases  undistorted.  The  range  of  observed  gains,  due  to  normal 
statistical  fluctuations  as  well  as  variations  in  background  noise,  is  shown  in  the  Table.  The 
range  for  delayed  summing  to  enhance  surface  waves  included  other  events  in  addition  to  th£ 
three  mentioned  above.  Figure  12  shows  the  waveforms  obtained  on  the  three  orthogonal  sets 
of  seismometers  in  processing  one  event  to  enhance  S  relative  to  noise.  Figure  13  shows  the 
results  of  processing  which  made  an  initially  invisible  LP  P-phase  quite  visible  on  the  processed 
traces.  Figure  14  shows  the  power  spectral  density  of  a  typical  vertical  seismogram  and  a  typical 
curve  of  SNR  gain  vs  frequency  for  maximum-likelihood  processing.  Curves  for  horizontal  com¬ 
ponents  are  similar. 

An  experiment  was  performed  to  determine  if  filters  which  are  designed  to  suppress  a  long- 
period  interfering  teleseism  on  one  particular  day  can  be  used  to  suppress  a  teleseism  from  the 
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TABLE  VI 

LONG-PERIOD  SIGNAL-TO-NOISE  GAINS  (db) 

(N  =  19,  v/N=  13  db) 

Type  of  Processing 

Component 

DS 

WDS 

FS 

P-phase 

Z 

8-11 

12-14 

13-16 

X 

8-13 

11-14 

12-16 

Y 

7-13 

11-14 

13-15 

S-phase 

Z 

8-11 

12-14 

13-16 

X 

8-13 

11-14 

12-16 

Y 

7-13 

11-14 

13-15 

Surface-wave  phase 

Z 

6-11 

Not  measured 

X 

10-12 

Not  measured 

Y 

9-11 

Not  measured 

same  region  on  another  day.  The  filter  coefficients  for  FS  and  WDS  were  obtained  by  placing  a 
21  November  1966  teleseism  from  the  Kurile  Islands  region  in  the  fitting  interval,  which  is  used 
to  design  the  filters.  Thus,  these  filters  were  designed  to  suppress  this  teleseism  and  to  pass  S 
waves  from  another  event.  The  filters  were  then  applied  to  a  12  November  1966  teleseism  from 
the  Kurile  Islands,  and  the  overall  suppression  obtained  by  FS  was  about  23  db,  although  the  sup¬ 
pression  at  times  rose  to  as  high  as  31  db. 

Program  limitations  in  the  frequency-domain  maximum-likelihood  program  did  not  allow  us 
to  apply  WDS  and  FS  processing  to  detect  surface  waves  using  the  full  LASA  aperture.  However, 
delayed  sum  processing  has  been  used  to  assist  in  the  analysis  of  long-period  surface  waves  ar¬ 
riving- at  an  array.  On  20  November  1966,  there  was  a  series  of  three  small  earthquakes  off  the 
Northern  California  coast  and  a  small  earthquake  (m^  ^4.5)  near  Rat  Island.  A  beam  steered 
toward  the  California  events  passed  these  surface  waves  while  rejecting  the  Rat  Island  surface 
waves  by  10  db.  When  the  beam  was  steered  toward  Rat  Island,  the  Northern  California  off¬ 
shore  event  surface  waves  were  attenuated  by  8db.  The  processing  thus  allowed  identification 
of  the  Rat  Island  surface  waves.  This  could  not  be  done  from  the  raw  traces.  The  beams,  as 
well  as  several  raw  traces,  are  shown  in  Fig.  15.  The  WDS  and  FS  suppression  of  surface  waves, 
described  above,  holds  promise  as  an  even  more  powerful  tool  than  delayed  summing  for  the 
identification  of  one  surface-wave  source  in  the  presence  of  another. 

J.  Capon 

R.  J.  Greenfield 

R.  T.  Lacoss 
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Fig.  12.  Three-component  LP-processed  traces  for  the  21  November  1966  Kurile  Islands  event 
(SP  P-magnitude  =  6.0,  LP  Rayleigh  magnitude  =  y 5 ). 
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Fig.  13.  Results  of  processing  the  12  November  1966  Kurile  Islands  event  to  uncover  hidden  LP  P-wave 
(SP  P-magnitude  =  5.1,  LP  Rayleigh  magnitude  =  5.3). 
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Fig.  14.  Typical  SNR  gain  vs  frequency  and  typical 
input  power  spectrum  for  vertical  long-period  array 
(21  November  1966  Kurile  Islands  event). 
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Fig.  15.  Beams  formed  with  long-period  vertical  seismometers  (20  November  1966  events). 
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V.  MONTANA  LASA  SYSTEM 


A.  AUTOMATIC  MONITORING  AND  FAULT  DIAGNOSIS 

In  the  previous  Semiannual  Technical  Summary  Report  (30  June  1966),  the  background  lead¬ 
ing  to  this  series  of  experiments  was  discussed.  A  brief  description  of  the  initial  goals  and  (he 
means  of  achieving  them  were  included.  System  details  were  not  available,  since  design  had  been 
under  way  for  only  a  few  months. 

At  present,  system  design  is  complete  and  hardware  installation  is  expected  to  be  complete 
within  a  month.  The  hardware  consists  of  an  SDS-93  0  general-purpose  digital  computer  coupled 
to  the  LASA  signal  collection  system  and  to  the  PDP-7  computer  by  means  of  a  Lincoln-designed 
interface  (MINS).  The  peripheral  equipments  initially  associated  with  the  SDS-930  are  largely 
commercial  products.  They  include  three  magnetic  tape  drives,  a  high-speed  line  printer,  £ 
card  reader,  a  paper-tape  reader/punch,  a  keyboard/printer  and  an  X-Y  plotter.  This  hardware 
is  currently  being  installed  at  the  LASA  Data  Center  at  Billings.  The  intent  in  pursuing  this  sys¬ 
tem  design  has  been  to  provide  as  much  flexibility  as  is  economically  feasible,  so  that  a  series 
of  experiments  can  be  performed  in  quick  succession.  Figure  16  shows  the  monitoring  system 
structure  in  block  diagram  form.  The  system  will  perform  two  distinct  functions: 

(1)  Momtor  Function  —  Checks  continuously  and  automatically  on  the  perform¬ 
ance  status  of  each  part  of  the  Montana  LASA,  reporting  to  the  PDP-7 
computer  and  to  the  maintenance  supervisor  when  performance  of  a  part 

is  subnormal. 

(2)  Diagnostic  Function  —  Enables  the  maintenance  supervisor  to  efficiently 
couple  with  the  Montana  LASA  system  for  purposes  of  making  diagnostic 
tests  and  measurements  on  system  elements  whose  performance  has  been 
questioned  by  the  monitor  function. 

The  initial  set  of  computer  programs  to  perform  these  functions  is  organized  as  an  executive 
program  controlling  the  operations  of  the  real-time  monitor  program  plus  the  diagnostic  pro¬ 
grams.  Figure  17  shows  the  timing  in  the  50-msec  data  cycle. 

The  executive  program  performs  all  input  and  output  operations,  monitors  and  services  all 
interrupts,  calls  in  and  executes  all  diagnostic  programs,  and  in  general  supervises  the  opera¬ 
tion  of  the  entire  programming  system.  Included  are  a  loader,  an  octal  dump  routine,  initializa¬ 
tion  routines,  a  recovery  routine  to  handle  power  failure,  a  routine  to  send  telemetry  comfnands, 
and  miscellaneous  subroutines  of  a  utility  nature. 

The  function  of  the  monitor  is  to  determine  how  well  the  Montana  LASA  hardware  system  is 
performing.  It  will  detect  and  identify  deficiencies  and/or  malfunctions  in  the  system  and  alert 
the  ooerator  and/or  the  PDP-7  to  the  particular  condition.  The  monitor  will  perform  the  fol¬ 
lowing  functions: 

(1)  For  each  sensor,  compute  the  absolute  average  of  sensor  data  over  some 
period  of  time  (minutes)  and  determine  if  the  average  exceeds  some  pre¬ 
viously  determined  threshold  value.  The  threshold  for  each  subarray  will 
be  determined  from  the  average  of  the  25-sensor  averages. 

(2)  For  each  data  circuit,  count  the  number  of  sync  word  errors  and  deter¬ 
mine  whether  the  allowable  error  rate  has  been  exceeded  over  some  unit 
time. 

(3)  For  each  data  circuit,  count  the  number  of  parity  errors  and  determine 
whether  the  allowable  error  rate  has  been  exceeded  over  some  unit  time. 
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Fig.  16.  Automatic  monitoring  and  fault  diagnosis  system. 
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Fig.  17.  Basic  data  cycle. 
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(4)  Maintain  subarray/sensor  status  for  transmission  to  the  PDP-7  when  a 
change  in  status  occurs  or  when  requested  by  the  operator. 

(5)  Periodically  sense  a  series  of  input  lines  to  check  the  on/off  status  of 
the  Data  Center  equipment,  and  a  non-comparison  of  the  two  Data  Center 
clocks. 

(6)  Monitor  the  21  telemetry  channels  for  alarm  conditions. 

The  output  of  the  monitor  will  take  several  forms.  Statistical  data  in  the  form  of  the  abso¬ 
lute  average  of  sensor  data  and  data  circuit  sync  and  parity  errors  will  be  output  to  the  paperi- 
tape  punch  every  two  hours.  Weather  information  obtained  from  site  AO  will  also  be  output  tq  the 
paper-tape  punch  periodically.  All  communications  between  the  program  and  the  operator  will  be 
reproduced  on  the  paper  tape  for  later  off-line  printing. 

Error  conditions,  system  malfunctions  and  deficiencies,  etc.,  detected  in  (1)  through  (6) 
above,  will  be  logged  on  the  console  typewriter.  All  such  printouts  will  be  accompanied  by  an 
audible  alarm. 

Several  diagnostic  programs  are  being  written.  The  programs  will  reside  on  the  system 
master  tape  and  will  be  available  to  the  operator  on  request  via  the  console  typewriter.  These 
programs  will  provide  the  operator  with  a  means  of  taking  a  closer  look  at  any  part  of  the  sys¬ 
tem  he  wishes  to  investigate.  The  diagnostic  function  will  be  performed  in  two  modes  —  manual 
and  automatic.  The  automatic  mode  will  be  set  when  the  system  is  loaded  and  can  be  reset  by  the 
type-in  "AUTO'1  at  any  time.  In  this  mode,  a  "resident"  diagnostic  will  be  in  core  and  the  mon¬ 
itor  can  perform  certain  diagnostic  functions  automatically.  The  operator  places  the  diagnostic 
function  in  the  manual  mode  by  calling  in  any  other  diagnostic  from  the  system  master.  Thd  di¬ 
agnostic  functions  presently  being  considered  are  those  which: 

(1)  Print  out  selected  data  compiled  by  the  monitor,  including  sensor  aver¬ 
ages  and  sync  and  parity  error  counts. 

(2)  Plot  sensor  data  on  the  X-Y  plotter. 

(3)  Transmit  and/or  receive  various  known  words. 

(4)  Exercise  the  test  SEM  at  Miles  City  (site  22)  to  verify  that  new,  modified, 
or  repaired  components  are  functioning  properly. 

(5)  Obtain  a  read-out  of  the  conditions  which  caused  an  alarm  indication  in  the 
telemetry  channel. 

(6)  Compute  and  output  the  DC  offset  of  the  multiplexer  at  each  site. 

(7)  Perform  a  Fourier  analysis  of  the  1-Hz  sine  wave  test  signal. 

(8)  Perform  a  Fourier  analysis  of  the  random  sequence  generator  test  signal 
to  obtain  the  transfer  function  of  a  selected  sensor  (broadband  calibration). 

J.  R.  Brown  E.W.  Richards 
J.  H.  Helfrich  R.  V.  Wood,  Jr* 

B,  LONG- PERIOD  SEISMOMETERS 

Installation  of  the  Montana  LASA  long-period  seismometers  was  completed  by  Geotech  Divi¬ 
sion,  Teledyne  Industries,  and  formally  accepted  by  the  Air  Force  Technical  Applications  Cen¬ 
ter  (AFTAC)  at  the  end  of  October  1966.  At  that  time  Lincoln  Laboratory  took  over  responsi¬ 
bility  for  the  long-period  system. 
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Fig.  18.  Long-period  system  transfer  function  (25-sec  filter). 
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Detailed  observations  of  long-period  data  have  been  conducted  only  for  a  month  or  two,  so 
that  relatively  little  information  on  component  performance  is  available.  The  intent  is  to  operate 
the  long-period  system  for  several  more  months  with  no  changes  in  parameters,  so  as  to  provide 
a  better  basis  for  performance  judgments.  With  this  in  mind,  the  2 5 -sec  filters  were  installed 
in  all  63  long-period  channels  during  October  1966  and  are  still  in  place.  The  nominal  system| 
amplitude  response  for  this  configuration  is  shown  in  Fig.  18.  Wideband  filters,  which  will  ex* 
tend  the  response  toward  1.0  Hz,  are  under  development  by  Texas  Instruments  and  will  probably 
be  installed  in  the  spring.  j  p  Densler 

Ft.  V.  Wood,  Jr. 


C.  SYSTEM  UPGRADING 

1.  Well-Head  Vault  Electronics  Package 

As  part  of  a  continuing  program  of  demonstrating  on  one  subarray  techniques  for  upgrading 
component  performance  in  the  entire  LASA,  a  completely  revised  well-head  electronics  pack¬ 
age  has  been  installed  in  subarray  E3.  Field  measurements  had  indicated  that  major  sources*  of 
cross-talk  and  coupling  between  signal  lines,  calibration  lines,  and  60-Hz  ground  currents  were 
inherent  in  the  original  well-head  vault  electronics  and  cabling.  The  new  installation  consist^  of 
two  separate  sealed  containers.  One  box  contains  all  cable  interconnections,  plus  lightning  pro¬ 
tection  circuits.  Normally,  entry  to  this  junction  box  after  installation  will  not  be  necessary. 
Another  box,  coupled  to  the  cable  junction  box  by  connectors,  contains  the  electronics  and  will  be 
replaced  as  a  unit  in  the  field  in  the  event  of  failure.  This  new  electronics  box,  as  installed  at 
E3,  includes  the  RA-5  redesign  to  reduce  gain  sensitivity  to  temperature.  An  important  feature 
of  the  new  well-head  vault  system  is  that  gain  adjustment  is  now  performed  at  the  central  ter¬ 
minal  vault  for  all  sensor  channels.  j  H  Helfrich 

C.  B.  Swanton 


2.  Extension  of  Subarray  E3 


Results  of  subarray  signal  processing  reported  in  Sec.IV-A  indicated  the  possibility  that 
signals  from  a  subarray  of  about  20-km  aperture  could  be  processed  on  line  with  a  signal-to- 
noise  improvement  considerably  closer  to  that  obtainable  from  off-line  maximum-likelihood  proc¬ 
essing  than  is  the  case  for  the  7 -km  subarray.  This  has  important  implications  for  future  LASA 
design.  To  verify  these  results,  and  to  test  engineering  feasibility  of  large  subarrays,  it  was 
decided  to  modify  subarray  E3  to  give  a  19-km  aperture,  with  25  sensors  and  a  minimum  sensor 
spacing  of  3  km.  This  modification  was  satisfactorily  completed  by  the  contractor  in  Montana 
in  early  December,  within  two  weeks  of  the  time  predicted  when  the  work  started  in  early 
September.  R.  Wood<  Jr. 


3.  Snbarray  SEM  Control/Output  Modifications 

For  six  or  eight  months,  it  has  been  clear  that  major  changes  in  the  SEM  telemetry  com¬ 
mand  sec:ion  were  required  to  adequately  meet  changing  needs.  Specifically,  a  major  source 
of  errors  in  the  digitized  data  arriving  at  Billings  was  found  to  be  the  improper  response  of  the 
SEM  telemetry  receiver  to  noise  on  the  Billings -to -subarray  channels.  Also,  the  capacity  of 
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the  telemetry  command  channel  had  been  reached,  and  could  not  readily  be  expanded  without  re¬ 
design.  The  eventual  control  of  the  long-period  seismometer  system  from  Billings  will  require 
a  substantial  number  of  new  telemetry  messages. 

It  was  decided  that  the  above  needs  could  best  be  met  by  designing  and  building  a  completely 
new  SEM  control  drawer,  and  extensively  modifying  the  output  drawer.  The  modification  provides 
a  telemetry  command  receiver  which  has  protection  against  errors  due  to  input  noise.  It  also 
provides  the  capability  for  63  distinct  telemetry  commands,  about  three  times  the  previous  ca¬ 
pability.  All  the  hardware  for  this  modification  has  been  completed  and  tested,  and  installation 
is  currently  under  way.  j  p  j}ensier 


36 


UNCLASSIFIED 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  classification  of  title,  body  of  abetract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 

1.  ORIGINATING  ACTIVITY  (Corporate  author) 

Lincoln  Laboratory,  M.I.T. 

2a.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2b.  GROUP 

None 

3.  REPORT  TITLE 

Semiannual  Technical  Summary  Report  to  the  Advanced  Research  Projects  Agency 
on  Seismic  Discrimination 

4.  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

Semiannual  Technical  Summary  Report  -  1  July  to  31  December  1966 

5.  AUTHOR(S)  (Last  name,  first  name,  initial) 

Green,  Paul  E. 

8.  REPORT  DATE 

31  December  1966 

la.  TOTAL  NO.  OF  PAGES 

44 

7b.  NO.  OF  REFS 

8 

8a.  CONTRACT  OR  GRANT  NO. 

AF  19(628)-5167 

b.  PROJECT  NO. 

ARPA  Order  512 

C. 

d. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S) 

Semiannual  Technical  Summary 
(31  December  1966) 

9b.  OTHER  REPORT  NO(S)  (Any  other  numbers  that  may  be 
assigned  this  report) 

ESD-TR-67-18 

10.  AVAILABILITY/LIMITATION  NOTICES 

Distribution  of  this  document  is  unlimited. 

11.  SUPPLEMENTARY  NOTES 

None 

12.  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects  Agenpy, 
Department  of  Defense 

13.  ABSTRACT 

Studies  have  continued  of  the  ability  of  a  single  Large  Aperture  Seismic  Array  (LASA) 
station  to  detect  and  roughly  locate  teleseismic  events  and  to  generate  outputs  derived  from 
various  long-  and  short-period  seismogram  parameters,  particularly  those  related  to  exact 
epicenter  location  and  source  discrimination.  The  behavior  of  networks  that  exchange  such 

1  data  and  that  consist  of  both  small  and  large  array  stations  has  begun  to  receive  attention. 

Computer-controlled  automatic  monitoring  and  maintenance  of  large  seismometer  arrays 
are  discussed. 

14.  KEY  WORDS 

seismic  array  seismometers 

seismology 

37  _ UNCLASSIFIED 

Security  Classification 


Printed  by 

United  States  Air  Force 
L.  G.  Hanscom  Field 
Bedford,  Massachusetts 


% 


